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Abstract 

Closed-form analytic solutions useful for the design of capillary flows in a variety of containers 
possessing interior comers were recently collected and reviewed. Low-g drop tower and aircraft 
experiments performed at NASA to date show excellent agreement between theory and experiment for 
perfectly wetting fluids. The analytical expressions are general in terms of contact angle, but do not 
account for variations in contact angle between the various surfaces within the system. Such conditions 
may be desirable for capillary containment or to compute the behavior of capillary corner flows in 
containers consisting of different materials with widely varying wetting characteristics. A simple 
coordinate rotation is employed to recast the governing system of equations for flows in containers with 
interior corners with differing contact angles on the faces of the corner. The result is that a large number 
of capillary driven corner flows may be predicted with only slightly modified geometric functions 
dependent on comer angle and the two (or more) contact angles of the system. A numerical solution is 
employed to verify the new problem formulation. The benchmarked computations support the use of the 
existing theoretical approach to geometries with variable wetability. Simple experiments to confirm the 
theoretical findings are recommended. Favorable agreement between such experiments and the present 
theory may argue well for the extension of the analytic results to predict fluid performance in future large 
length scale capillary fluid systems for spacecraft as well as for small scale capillary systems on Earth. 


Introduction 

A variety of closed form solutions are available that predict important features of capillary driven 
flows along interior comers (ref. 1). Such flows are common in porous media and micro-fluidic devices, 
but are especially prevalent, and over large length scales, in the fluids management systems of spacecraft 
(ref. 2). For example, concerning the latter, interior corners are a common passive geometric solution to 
the positioning of liquid fuels and cryogens in their respective containers as shown in figure 1. In this 
figure the interior comers formed by vanes within a “spheriodal” tank (refs. 3 and 4) are shown to 
maintain the fuel over the outlet to the engine for proper operation during in-flight maneuvers in a low 
gravity environment. 

Important characteristics of such flows include interface configurations, stability, and flow rate with 
respect to time and system perturbation. Solutions for flows in cylindrical polygonal geometries (refs. 5 
to 8), planar and non-planar ‘infinite’ comers (refs. 9 to 1 1) and three-dimensional geometries (refs. 12 
and 13) have been reported. 

In this paper the generalized approach of previous analyses (ref. 8) is extended to capillary driven 
flows within planar interior comers where the wetting properties of the fluid differ between the two faces 
that form the interior comer. Such flows arise in systems consisting of differing materials, coatings, or 
other surface treatment. Several example problems are sketched in figure 2. 
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Figure 1 . — Heavily vaned VTRE tank (refs. 3 and 4) schematic showing 
fluid location in low-g environment. Ta nk is partially liquid 
filled as shown in a. with fluid occupying shaded region. 



Figure 2. — Example container sections with variable wetabilities on various surfaces. 


The differing wetabilities on each surface of the interior comer are characterized by contact angles 9i 
and 0 2 . This work is pursued with the hope of providing simple geometrically dependent coefficients that 
can be used to modify the closed form solutions of previous investigators for a variety of problems as 
cited above. 


Review of Problem Formulation 

A schematic of a typical interior comer flow problem is sketched in figure 3 with characteristic 
dimensions noted for a system where 9i = 0 2 . This is the case of previous investigations where the wetting 
characteristics are uniform between the fluid and walls of the comer producing a symmetric interface in 
the jc-v- plane (fig. 3b). 

In brief, by noting that the column of fluid in the comer is slender, namely ( H/L ) 2 = s 2 c 1 , and that 
surface curvature (~ 1 IHf) in the x-y-plane dominates the capillary pressure (s 2 / « 1), where /is the 
measure of interface curvature, the free surface S(y,z,t ) may be expressed as a construct of circular arcs in 
x-y-planes that satisfy the slope condition (i.e., contact angle condition) at the wall. The Navier-Stokes 
equation reduces to the z-component equation only, which when integrated numerically can be used to 
determine the z-component velocity distribution w(x,y,z,t). Numerically integrating this velocity over the 
cross-flow section (fig. 3b) yields the average velocity (w) = (w)(z,t ) , which when substituted into the 
global mass balance equation in the z-direction yields a governing partial differential equation for a 
variety of flows differing only by respective boundary conditions. Proper scaling of the governing 
equations significantly reduces the reliance on numerical data. 
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Figure 3. — Schematic of corner flow: a. characteristic dimensions, 
b. section detail, c. dimensionless section. 


Choosing scales x ~ H,y ~ H tana, z~L,S~H, P~o IHf, t ~ L/W, and 


s a sin 2 a 

w = W , 

h/ 


where 


( 1 ) 


/ = 


sma 


cos 9 -sin a 

The dimensionless problem outlined above for 8 2 c 1 and 8 2 /<C 1 yields free surface 

S(y,h) = h( 1 + f)-{f 2 h 2 -y 2 tan 2 a) I/2 , 
where y < \ f /;sin (5 /tan a| , z-component momentum equation 


• 2 2 
P 2 = h' a:x sin" a + vv lT cos a 


subject to no slip on the walls 


w - 0 on y = ±x , 

and no shear on the free surface 

2 

vv x - S v w v cot a = 0 on x = S , 

and the governing partial differential equation from global z-component mass balance, noting that 

P- - h-/h ? . 


( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
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(7) 


h t =F i {h;+hh zz ). 

The dimensionless cross-flow problem, equations (4) to (6), is sketched in figure 3c. The 
dimensionless geometric flow resistance coefficient is F, = F, (0,a) must be determined numerically. 
However, due to the choice of scales for this problem, in particular W, 1/8 < F, < 1/6 for all values of 0 
and a satisfying the critical wetting condition studied mathematically by Concus and Finn (ref. 14), 
where 0 < nl 2 - a in order to achieve the capillary flow condition studied here. For the remainder of this 
paper this condition shall be known as the C-F condition. If 0 > n!2 - a the liquid will not spread along 
the comer, but will retract until equilibrium is established. For many practical problems of interest F, may 
be treated as a constant, say F, ~ 1/7, and absorbed into the timescale of equation (7). The scale analysis 
method to determine W is presented in (ref. 10) and is instructive for subsequent equation 
nondimensionalization prior to asymptotic analyses. (The simple problem of fully developed laminar flow 
through a rectangular duct is revisited herein (Numerical Solution to Cross Flow Problem section) as a 
numerical benchmark and secondly to illustrate the utility of the general scaling approach. The hydraulic 
diameter and scale analysis solutions for average velocity are compared, the latter of which more 
narrowly brackets the exact solution for this problem.) The results of numerical solutions for F, (a,0) also 
demonstrate the value of a correctly scaled geometry. 


Formulation for Problem 0 t ^ 0 2 


A sketch of the dimensional comer flow cross-section for the two contact angle case is provided in 
figure 4a for a comer half-angle a. The comer is positioned symmetrically about the x’-axis, similarly as 
the comer in figure 3b. The C-F condition for this case is 


01+07 tc 

— — < a . 

2 2 


( 8 ) 


The further constraint of a partially wetting system, 0j < n/2 and 0 2 <n/2, places the interface 
behavior of interest here within the D\ -domain identified by Concus and Finn (ref. 15). The 
configuration of figure 4a is rotated in figure 4b such that the minimum height of the meniscus occurs 
along y = 0, which is to say that the radius of interface curvature in the x-y- plane pivots about a point on 
the x-axis whereas in figure 4a it does not. Thus, noting that f\ =f 2 =/ and oq + 0,2 = 2a, ai and a 2 may be 
computed from 


sin 2a 

T\ = tan a . = 

(N + cos 2a) 


and 


sin 2a 

T 2 =tana 2 = — — , 

N~ l + cos 2a 


(9) 


( 10 ) 


where 
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( 11 ) 


N = 


COS0 2 

COS0J 


To avoid redundancy in the analysis and computations to follow 0| < 0 2 and a.| > a 2 are defined. 
Thus, 0 < A < 1 and 0 < 0j / 0 2 < 1 . Note also that Sj = nil - - 0 t and 8 2 = n/2 - a 2 - 0 2 . By choosing 

alternative scales for y and w, namely; y ~ TH and 


w~W = 


so 

fi/ 


vl + T‘ y 


( 12 ) 


where 


T^Ti+Ti), (13) 

a convenient dimensionless cross flow problem is posed as sketched in figure 4c. In this figure the flow 
domain of figure 4b is stretched such that the dimensionless cross flow area is ~ (9(1) for all values of a, 
0i, and 02. The numeric problem of flow through the dimensional section sketched in figure 4b will be 
described shortly. 





Figure 4. — Schematic representations of cross flow sections with differing wetabilities on faces of interior 
comer, 0 b 0 2 : a. Comer symmetric about dimensional x’-axis, b. Axes rotated such that 
radius of interface curvature aligns with dimensional x-axis, c. Dimensionless, rotated section. 
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The dimensionless geometric scaling function T= T( a,0i,02) is used to account for the widely varying 
characteristic v-length scale arising for wide variations in a, 0i, and 02. Some limiting cases for T are 
provided below: 

1. Symmetry Condition, 0i = 0 2 . For this condition ai = a 2 = a and N = 1. Thus, from equation (13), 

2. T= tan a and the governing system given by equations (1) to (7) is recovered. 

3. Large Contact Angle, 0 2 — >n/2 with a < n/ 4 . For this condition, in order to satisfy the C-F 
condition (eq. (8)), 0i <c 1. Thus, N — * 0, a 2 — > 0, and a! — * 2a, and from equation (13), 

4. T— * (1/2) tan 2a. 

5. Small Comer Limit, a C 1. In this limit T~ a, which recovers the result of the symmetric scaling 
resulting in equation (4). 

6. Large Comer Limit, a 2 — » n/2. Defining Q = nil - a, this limit is identical to Q -< 1. Equation (13) 
reduces to T~ 1/Q when fi< 1, which also recovers the result of the symmetric scaling resulting in 
equations (4) to (6). 

Thus, the limiting behavior of T for the 0i ^ 0 2 problem suggests that F, computed numerically should 
be a narrowly bound function just as when computed using the symmetric scaling quantities for the 0i = 
0 2 problem. 


Modified System for Problem 0! ^ 0 2 

The governing dimensionless system of equations again constrained by s 2 <C 1 and s 2 /<C 1 yields 


/i- 


sin tt| 

COS0J - sin a! 


= / = 


sin a 2 

cos0 2 - sin a 2 


= f 2 ■ 


free surface 


S(y,h) = h(l+f)-(f 2 h 2 -T 2 y 2 ) 112 , 


where sin 8j < Ty / fh < sin 5 2 , z-component momentum equation 


P z = 


( j2 ^ 

\ l + T ‘ J 


+ 


' 1 ^ 
vi+r j 


w 


yy- 


subject to no slip on the walls 


„ xT, 

w = 0 on v- 

T 


w = 0 on y = ~ 


xT~, 


(14) 


(15) 


(16) 


(17) 

(18) 


and no shear on the free surface 
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on 


x = S , 


(19) 


w x~ T 2 S y w y =0 

and again the governing partial differential equation from global z- component mass balance is 
equation (7). 


Analytic Solutions to Cross-Flow Problem 

The ‘cross-flow’ problem described by equations (16) to (19) may be solved analytically in the limits 
of large and small interior comer angle. Introducing Q = n/2 - a, these limits are described 
asymptotically by Q 2 1 and a 2 C 1, respectively. Using normal series expansions of the dependent 
variables ( w , P, and S) in powers of either Q 2 or a 2 , leading order expressions for F, may be determined 
for comparisons with subsequent numerical calculations. Some important details of the asymptotic 
solution approach are provided here. 


Large Corner Angle Solution, Q 2 <sc 1 

The large angle limit is the condition a — > til 2 , or preferably Q — > 0 . Choosing expansions 


2 

w-w 0 + W] + ... 

(20) 

P = P 0 +Q 2 P x + ... 

(21) 

S - S 0 +Q 2 5j +..., 

(22) 


for conditions where Q 2 <c 1, the governing z-component momentum equation (16) reduces to 

Pz =w xx + 0(Q 2 ). (23) 

The no slip boundary conditions on the walls, equations (17) and (18), remain unchanged, but the free 
surface shear stress condition equation (19) simplifies to 

w x = 0+O(Q 2 ) on x-S 0 (24) 

Because P : = hJh 1 is not a function of v the velocity profile may be determined analytically from (23) 
producing 




(1 + T 2 )h z 


T 2 h 2 


PP__ S 

^2 ° T* 


2 P 


! i 


(1 + T 2 )h z + 

—w 

T 2 h 2 


( 25 ) 


for 0 < y < y m i and 
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(26) 


w 


°2 


(1 + T l )h z 


T 2 h 2 


T -X-s„i- T -^ 

2 T? T 2 


(\ + T 2 )h z + 

~ — ^ — w 

T 2 h 2 ° 2 


for y m2 < y < 0 . The two flow regions are depicted in figure 4c. The average velocity over the cross-flow 
section may be determined by integrating equations (25) and (26) over their respective sections such that 


(■ W o) = 


(1 + T-)h z 


Th 4 F A 


C f w + dxdy+ [ f w + dxdv = —F i h z 
Jo STy!T x °\ " J v„ i2 J-7)’/T ° 2 


(27) 


where y mi = (-l) 1+! h f sin(Q, - 0 ; )/7' . The dimensionless area function F A is given by 


F a ={F a +F Ai )I 2, 


(28) 


where, 


F Ai =r 


cos0 ; - sin(I2 ; - 0j-) 
cosQ, 


-(n f -0 f ) 


(29) 


The analytic solution to equation (27) is simplified significantly by first writing each £2, and 0, in terms of 
Q and then expanding all geometric functions for small Cl. The following definitions are employed to this 
effect: 


= 2 m £2/(1 + m), £2 2 = 2 £2/(1 + u) , where u = £2j /£2 2 with £2j < £2 2 , (30) 

0j = (1 + m) 0j / 2 £2 , 0 2 = (1 + m)Q l / 2 £2 , where m = 0j /0 2 , 0j < 0 2 , (31) 


and thus 


0 t = 2 m 0 2 £2/(1 + m), 0 2 = 20 2 £2/(1 + m) . 

sfc 

Noting that f=f u may be determined analytically as a function of 0 2 and m; 


(32) 


u — 


1 — (1 — m 2 ) 0 2 

1 + (1 - » 7 2 ) 0 2 


+ 0(Q 2 ) . 


(33) 


With knowledge of u, equations (25) and (26) may be substituted into the terms present in equation (27) 
and the integration performed to reveal 


Fj - fs + 20ni\ 0 2 + 8 (l + 25 m / 2 + nF )o 2 ~ — 6 8 m 1 (l - 50 m / 17 + m‘ 


)ef 


- 70 m 2 (1 + Am + m 2 ) 0 


*4 

2 
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+ 140777! 777 9 02 + 140 777 777 2 0 2 — 2 8 777 j 777 o 0 2 ~ 7 777 2 777 f 0 2 j 
^35^1 - 777 f0 2 2 j (l + 777 2 0 2 )^l + 2777 1 0 2 -3777?0 2 2 j + (9(Q 2 ), 

where for brevity in this expression ni \ =1 + 772 and 777 2 = 1 - m . Thus, under the constraint Q 2 <c 1 , 

* 

Fj — Fj (777, 0 2 ) . Some important limiting values of F, are: 


1. Symmetric wetting condition (9i = 0 2 ): m = 1 


^•( 1 , 02 ) 


5 + 4(10 + 290 2 +320 2 2 )0 2 
35(1 + 20 2 ) 2 (1 + 40 2 ) 


+ 0(C1 2 ), 


which agrees with previous analysis for this condition (ref. 16). Note that F, (1,0) = 1/7. 


(35) 


2. The C-F condition (Oj + 0 2 = 2 ST, or 0 t + 0 2 + 2a = 7r ): 0 2 =1/(1 + 777 ) 

F t =i + 0(fi 2 ). (36) 

6 

For this flat interface configuration F, is identical to the value determined analytically for the symmetrical 
problem for the C-F condition (ref. 16). 


1 . Perfect wetting ( 0 1 = 0): m = 0 


F,W 2) J + 3 ^ +9 < + 498» 2 4 +0(fl2) 

35(1 + 0 2 ) 4 (1 + 30 2 ) 


(37) 


For 0 2 = 1,F,(0,1)= 1/7. Thus for all combinations of 777 and 0 2 , 1/7 < F ; (t? 7,0 2 ) < 1/6 whence 1. 


Small Corner Angle Solution, a 2 C 1 

A similar asymptotic solution approach may be used to compute F t for small comer angles where 
a 2 -< 1. Under these constraints equation (16) reduces to P : = w vv +0(a 2 ) and applying the appropriate 
boundary conditions leads to 




(1 + T 2 )h z 
2T 2 h 2 



+ T(T l -T 1 )xy 



which applies for all y in the domain. Integration over the section yields 


( 38 ) 
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( 39 ) 


Fj — h Ot A(0j , 0? ) + 0{d ) 

6 

where F(Q l , 0 2 ) ~ 0(1) and is a known but cumbersome geometric quantity. It is sufficient here only to 
present the limiting values of F,- determined in this manner, and from equation (39), F t =1/6 + 0(a)for all 
values of 0i and 0 2 . 


Numerical Solution to Cross Flow Problem 

The analytic results for F, in the extreme limits of a, 0i, and 0 2 offer compelling evidence that F, is 
again a narrowly bounded function ~ 0(1/7) < F, < 1/6 for the 0i ^ 0 2 problem. However, for the majority 
of corner flows that do not appear to yield analytic solutions, equation (16) is solved numerically. The 
investigation is undertaken to verify the scaling employed for the 0i ^ 0 2 problem leading to the system of 
equations (14) through (19). If F, can be confirmed to be bounded for all a, 0i, and 0 2 , similarly as for the 
case 0i = 0 2 , then previous analytic solutions for a variety of problems can be resolved simply by 
replacing the flow resistance coefficient F, sin 2 a with F, T 2 /( 1 + T 2 ). 

A Mat lab’ code written for this purpose is first benchmarked against the geometrically simpler 
problem of fully developed laminar flow in a rectangular duct for which an analytic solution is available 
for comparison. This problem serves to illustrate both the subtlety and importance of the non- 
dimensionalization for the numerical approach. 

The system of equations (16) to (19) is then solved numerically and benchmarked against limiting 
values for the 0i = 0 2 problem that are analytic and/or have been computed by previous investigators 
(refs. 17 and 18). The term selected for comparison to the analytic results is the dimensionless flow 
resistance coefficient F, which for the dimensionless problem provides proportionality between average z- 
component velocity and interface slope, or 


(w) = -F i h z • (40) 

[Note: The viscous resistance of the flow is actually inversely proportional to the nearly constant function 
A,.] The numerically computed values for (w) are converted to F, for the comparisons. 


Benchmark: Fully Developed Laminar Duct Flow 

The classic problem of single -phase (no interface) fully developed laminar flow through a rectangular 
duct of half width a and half height b is selected for use as a benchmark for the code. An exact solution to 
Poison’s equation for this flow is available (ref. 19). The analytic equation for the dimensional average 
velocity is 


(w) 

\ t exact 


APa 2 
3 pL 


1- 


192 


I 


tanh(mq/2) 


K h ;=1,3,5,... 


• 5 


(41) 


where, q = b/a. Choosing length scales x~ a,y~ b an identical scaling procedure as employed in the 
development of governing equations (14) and (16) produces the velocity scale, 


W IFf-ect 


( A 


(1 + q 2 ) 


(42) 
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where, P z = AP/ L . Using this velocity scale, Poison’s equation is nondimensionalized yielding the 
governing equation, 


1 = 


( 7 A 


( ! ) 

n" 

w xx + 

1 



li+q “ ) 


l 1+ q J 




( 43 ) 


The modified Poison equation (43) for w(x,y) is solved numerically and integrated over the domain to 
determine the nondimensional average velocity, which is redimensionalized for comparisons to the exact 
solution of equation (41). The computational error for representative values of q is displayed in figure 5. 
As observed from the figure, acceptable convergence is achieved when 50,000 or more elements are used. 
Contour plots for q = 0, q = 1, and q — > <x> are presented in figure 6. The scaling permits accurate 
numerical solutions in the hard limits of q = 0 and q — > go. The value of the present scaling ( W rect , 
eq. (42)) is further illustrated by comparing to solutions obtained employing a hydraulic diameter scaling 
where x ~y ~ D hyd = abt(a + b), defined such that limiting values of are equivalent for comparisons in 
figure 7), 


^ hyd 


AP b 2 

(1 + q ) 2 


(44) 


4.0 

3.5 

g 3.0 

§ 2.5 
W 

§ 2.0 
.o 

3 1.5 


o 

U 


1.0 

0.5 

0.0 
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0 

□ 

gf 

A 

□ 

S8 

) A 
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[# 

P ^ 

A 

o □ 
CD 


□ 

o 


n 

0 5.45 
□ 96.8 
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x 6974 
o 12347 






0 20 40 60 80 100 120 


# of Elements (TO 3 ) 

Figure 5. — Convergence study of fully developed laminar flow through a 
rectangular duct. The relative error calculation compares equation (41) 
to numerical solutions calculated with the Matlab® PDE 
Toolbox for various selections of q. 
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Figure 6. — Velocity contour maps in the scaled domain for rectangular duct flow. 
The left profile is for r| = 0, the middle is q = 1 , and the right is q — » oo. 



n 


Figure 7. — Comparison of exact solution scaled by equation (42) 
and (44) respectively. The variation of <w> 
for the former is significantly reduced. 


The exact analytical solution for the average velocity through the rectangular duct equation (41) is 
normalized by both equation (42) and (44) scaling and presented in figure 7. The figure displays the 
extent to which the hydraulic diameter scaling increases the reliance on ‘numerical data’ over the full 
range of q, namely; 1/3 < {w) hyd < 9/16, whereas 0.281... < {w) rect < 1/3. Thus, W rect captures a higher 
degree of the flow’s dependence on duct geometry and confines solutions for all q to a domain over 
40 percent smaller than for the otherwise acceptable hydraulic diameter scaling. The scaling method for 
W, equation (12), for the comer flow problem addressed herein plays a similar role. [Note: Other methods 
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to nondimensionalize this simple rectangular duct flow problem can lead to situations that are 
significantly worse than even the hydraulic diameter scaling. For example, Patzek and Silin (ref. 20) 
employ a dimensionless flow conductance g where W PS ~ AP ab/\i L. When used to nondimensionalize 
equation (41), such a scaling does not produce results that are bounded by (9(1) constants, but rather 
0 < (w) < (9(1) for all h, 0 < p < oo.] 


Benchmark: Corner Flow with 0] = 0 2 

The Matlatf numerical procedure is benchmarked for the comer flow problem for several analytic 
solutions of Fj with ©j = 0 2 = 0. The favorable comparisons are shown in table 1 . For case III, a = 89.8° is 
used to compute the numerical value of F); for all other cases the computations are carried out exactly as 
stated. As was the case for the rectangular duct, the fact that numerical data can be computed in the 
extreme limits of a and 0, is a result of the scaling of the problem, which produces ~ (9(1) area domains 
for all values of0<a<7t/2, O<0 1 <7r/2 and 0 2 < n - 2a - ©j satisfying the C-F condition. 

When possible, and for the case 0 t = 0 2 , the numerical calculations are found to be in close agreement 
with the previous results of Ayyaswamy et al. (ref. 17) and Ransohoff and Radke (ref. 18). 


Table 1. — Comparison of analytical and numerical values of F t . 


Case 

Condition 

Equation 

F, 

analytic 

Ft 

numerical 

Error 

(%) 

I 

a =45° 
0 = 45° 

eq. (41) 
w/ q = 1 

0.140577 

0.140571 

0.0044 

II 

a = 0 

eq. (39) 

1/6 

0.166664 

0.0016 

III 

a — > 90° 
0 = 0 

eq. (35) 

0.142857 

0.142844 

0.0091 

IV 

a = 90° 
5 = 0 

eq. (36) 

1/6 

0.166651 

0.0094 


Numerical Results for F,- (a,6i,6 2 ) 

F,(a,0i,0 2 ) is determined numerically using Mat I ah" with PDE tool box in Microsoft Windows' 
on a 1 GHz machine with 512 Mb RAM. An adaptive mesh is employed containing a minimum of 
50,000 elements and requiring 5 to 15 minutes to converge for a given choice of a, ©i, 0 2 . figure 8 
presents over 2000 values of F, computed for all wetting conditions satisfying the C-F condition and 
where O<OJ<O 2 <TCZ 2 and 0<a<7c/2.1t should be noted that certain wetting conditions have not 
been considered in this paper; specifically 0 2 > 7c/2 with 0i < Tt - 2a - 0 2 . From figure 8 it is found that 
1/9 < Fi < 1/6 for all values chosen for a, 0i and 0 2 . As observed in Table 1, the computational errors are 
extremely low, < 0.5 percent for these calculations. All computed values for F,(a,0i,0 2 ) are tabulated in 
the Appendix. For the problem 0i = 0 2 , 1/8 < F) < 1/6. Thus, the scaling for the problem 0i ^ 0 2 problem, 
which recovers the results of problem 0i = 0 2 for symmetric wetting, yields a similar though slightly 
expanded range of values for F,. Without a priori knowledge of the specific value of F, one could assume 
Ft ~ 1/7.2 and be assured of at worst +20 percent accuracy. 
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Figure 8. — Scatter plot showing all the numerically calculated values 
of F’,(a,9i,02), where (9i + 02)/2 < tt:/2 - a, 9i < nil, 

02<7r/2, and0<0 1 /0 2 < 1. 




Figure 9. — Velocity contours in the scaled domain. The profile on the 
left is a = 0 and on the right is a = n/2. 



Figure 10. — Velocity contours in the scaled domain for a = n/4. The profile on the left is for the wetting 
condition 0! = 0 and 02 = 0, the middle is 0i = 0 and 02 = 7t/4, and the right is 0i = 0 and 02 = tt/2. 
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Application of Results 

The recent review of comer flows (ref. 1) provides a collection of closed form solutions for a fixed 
contact angle (9i = 0 2 ). The solutions are in part described by functions G and F A which are modified by 
the results of this work for the two contact angle problem: 


G = 


c ■ 2 
a Fj sm a 


F / 

F / 


\ + T 


(45) 


F A => F A 


= \( F Ai+ F A 2 ) 


(46) 


The dimensionless geometric viscous resistance function T 2 /( 1 + T 2 ) for the general 0! ^ 0 2 replaces 
sin 2 a for the 9i = 0 2 problem. Recall that T 2 /( 1 + T 2 ) = sin 2 a when 0! = 0 2 , since T= (tanoq + tana 2 )/2 
with ai and a 2 given in equations (9) and (10). The dimensionless cross flow area function F A represents 
the average of F A values computed using 0| and 0 2 . Note also that F A = F A when 0i = 0 2 . 

The more general forms for equations (45) and (46) may be simply substituted wherever they appear in 
(ref. 1) to determine the flow in comers for cases where 9j ^ 0 2 . 

From the values of F l (a, 9i, 0 2 ) listed in the Appendix it can be shown that errors in F, of less than 
5 percent are incurred using the simplifying approximation F, (a, 9i, 0 2 ) ~ F, (a, Q ave ) with Q ave = (9i + 
9 2 )/2. The conditions under which this approximation holds are 0 2 - 0i ~ n/6 and/or a + 0 2 < n/2, where 
0 < 0i < 0 2 < 7t/2 and 0 < a < nil. 

The special case of spontaneous capillary rise requires the calculation of a constant height Fl 
boundary condition in each comer of the container where the C-F condition is satisfied. This may be 
accomplished by extending the method of de Lazzer et al. (ref. 21) to compute the radius of curvature of 
the interface R: 


where 


R = fiHt = 


2X 


f 

1 - 

V 


4 AX 


\l/2 



(47) 


n 


P=^Pi COS0; , 
i= 1 

(48) 

2 = t F An , f , 
1=1 

(49) 

FAni=F Ai lff. 

(50) 


X is the total cross flow area function for the container and P is a section perimeter function weighted by 
cos0, for each face of width P,. A somewhat generalized polynomial container section for this flow 
scenario is represented schematically in figure 1 1 identifying the notation used in equation (47). When the 
C-F condition is not met in a particular comer i, F Ani = 0 for that comer. Using this approach a global 
similarity solution for flows in all comers satisfying the C-F condition is also possible when 0i ^ 0 2 
(refs. 4 and 7). 
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Figure 11. — n-sided polygonal section. 


Recommended Experiments 

Perhaps the simplest experiments capable of testing the theoretical work reported here are capillary 
rise drop tower experiments conducted in an identical manner as those performed by Weislogel and 
Lichter (ref. 8). In these experiments a right cylindrical container with interior comers is partially filled 
with a wetting liquid and place vertically in the drop tower test apparatus. Release of the experiment into 
free fall results in the spontaneous rise of liquid in the comers which is photographed at large working 
distances to accurately determine the interface location as a function of time. 

The cross section of the container is critical to gather highly quantitative data from the images. For 
example, a right equilateral triangular cylindrical container can be constructed for the two contact angle 
condition with a priori knowledge of the contact angles of the two surfaces. A ray trace analysis can then 
be performed such that mismatched refractive indices between container and test fluid are minimized and 
the meniscus height h(z,t) may be accurately imaged for comparison with predictions. In general, the 
containers should be designed such that the interface may be observed from the profile point of view of 
figure 4b rather than figure 4a. A sample sectional view of such a container is shown in figure 12. 
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Figure 12. — Equilateral triangular container section 
with differing contact angles for drop tower tests. 
Line of sight shown at right (not to scale). 


Concluding Remarks 

The interior comer flow problem (for both 9i = 0 2 and Cf ^ 0 2 ) appears to be correctly scaled by the y 
length scale y ~ TH. This scaling produces a dimensionless flow resistance coefficient F, that is calculated 
and computed herein and shown to be narrowly confined such that 1/9 < 1/6. These results bode well 

for the many previous comer flow solutions that need only replace a coefficient to extend their relevance 
to flows where the wetting conditions on opposing faces of the comer differ. 
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Appendix A — Tables 2 through 20 
Computed Values of F,(a,0i,0 2 ). 


The format for Tables 2 through 20 follows the form 


a 

02 


0i 

F,(a,0 h 0 2 ) 






Table 2. — -F,(a, 0 |, 02 ) for a = 0. 
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0 
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0.167 

0.167 
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0.167 

0.167 

0.167 

0.167 

0.167 
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5 
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10 
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0.167 
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15 
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0.167 
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0.167 
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0.167 
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0.167 
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0.167 
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35 
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0.167 
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40 
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0.167 
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45 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 

0.167 
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0.167 
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65 
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0.167 

0.167 

0.167 
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70 
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0.167 
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0.167 
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0.167 
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75 
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0.167 
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0.167 
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Table 3. — -F^a, 61 , 62 ) for a = 5. 
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Table 4. — -F,{a,0 1 ,0 2 ) for a = 10. 
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0.157 

0.157 

0.158 

0.159 

0.160 

0.161 

0.162 

40 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.155 

0.155 

0.155 

0.155 

0.156 

0.156 

0.157 

0.157 

0.158 

0.159 

0.160 

0.161 

0.161 

45 

0.154 

0.154 

0.154 

0.154 

0.155 

0.155 

0.155 

0.155 

0.155 

0.156 

0.156 

0.156 

0.157 

0.157 

0.158 

0.159 

0.159 

0.160 

0.160 

50 

0.155 

0.155 

0.155 

0.155 

0.155 

0.155 

0.155 

0.156 

0.156 

0.156 

0.156 

0.156 

0.157 

0.157 

0.158 

0.159 

0.159 

0.160 

0.160 

55 

0.155 

0.155 

0.155 

0.155 

0.156 

0.156 

0.156 

0.156 

0.156 

0.156 

0.156 

0.157 

0.157 

0.158 

0.158 

0.159 

0.159 

0.159 

0.159 

60 

0.156 

0.156 

0.156 

0.156 

0.156 

0.156 

0.156 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.158 

0.158 

0.159 

0.159 

0.159 

0.158 

65 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.157 

0.158 

0.158 

0.158 

0.159 

0.159 

0.159 

0.159 

0.157 

70 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.158 

0.159 

0.159 

0.159 

0.160 

0.159 

0.156 

75 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.159 

0.160 

0.160 

0.160 


80 

0.160 

0.160 

0.160 

0.160 

0.160 

0.160 

0.160 

0.160 

0.160 

0.159 

0.159 

0.159 

0.159 

0.159 

0.160 

0.160 

0.161 



85 

0.161 

0.161 

0.162 

0.162 

0.162 

0.162 

0.161 

0.161 

0.161 

0.160 

0.160 

0.159 

0.159 

0.159 

0.159 

0.160 




90 

0.163 

0.163 

0.163 

0.163 

0.163 

0.163 

0.162 

0.162 

0.161 

0.160 

0.160 

0.159 

0.158 

0.157 

0.156 






Table 5. — ■F,{a,0 1 ,02) for a = 15. 


15 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

0 

0.146 

0.146 

0.146 

0.146 

0.147 

0.147 

0.147 

0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.152 

0.153 

0.153 

0.154 

0.154 

5 

0.146 

0.146 

0.146 

0.147 

0.147 

0.147 

0.147 

0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.152 

0.153 

0.153 

0.154 

0.154 

10 

0.146 

0.146 

0.146 

0.147 

0.147 

0.147 

0.148 

0.148 

0.148 

0.149 

0.149 

0.150 

0.151 

0.151 

0.152 

0.153 

0.154 

0.154 

0.154 

15 

0.146 

0.147 

0.147 

0.147 

0.147 

0.148 

0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.152 

0.152 

0.153 

0.154 

0.154 

0.154 

20 

0.147 

0.147 

0.147 

0.147 

0.147 

0.148 

0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.152 

0.153 

0.153 

0.154 

0.154 

0.154 

25 

0.147 

0.147 

0.147 

0.148 

0.148 

0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.151 

0.152 

0.153 

0.153 

0.154 

0.154 

0.154 

30 

0.147 

0.147 

0.148 

0.148 

0.148 

0.148 

0.149 

0.149 

0.149 

0.150 

0.150 

0.151 

0.151 

0.152 

0.153 

0.153 

0.154 

0.154 

0.153 

35 

0.148 

0.148 

0.148 

0.148 

0.148 

0.149 

0.149 

0.149 

0.150 

0.150 

0.151 

0.151 

0.152 

0.152 

0.153 

0.153 

0.154 

0.154 

0.152 

40 

0.148 

0.148 

0.148 

0.149 

0.149 

0.149 

0.149 

0.150 

0.150 

0.150 

0.151 

0.151 

0.152 

0.152 

0.153 

0.153 

0.154 

0.153 

0.151 

45 

0.149 

0.149 

0.149 

0.149 

0.149 

0.150 

0.150 

0.150 

0.150 

0.151 

0.151 

0.152 

0.152 

0.153 

0.153 

0.153 

0.153 

0.153 

0.150 

50 

0.149 

0.149 

0.149 

0.150 

0.150 

0.150 

0.150 

0.151 

0.151 

0.151 

0.151 

0.152 

0.152 

0.153 

0.153 

0.153 

0.153 

0.152 

0.149 

55 

0.150 

0.150 

0.150 

0.150 

0.150 

0.151 

0.151 

0.151 

0.151 

0.152 

0.152 

0.152 

0.153 

0.153 

0.154 

0.154 

0.153 

0.152 

0.148 

60 

0.150 

0.150 

0.151 

0.151 

0.151 

0.151 

0.151 

0.152 

0.152 

0.152 

0.152 

0.153 

0.153 

0.154 

0.154 

0.154 

0.154 

0.152 

0.147 

65 

0.151 

0.151 

0.151 

0.152 

0.152 

0.152 

0.152 

0.152 

0.152 

0.153 

0.153 

0.153 

0.154 

0.154 

0.154 

0.155 

0.154 

0.152 


70 

0.152 

0.152 

0.152 

0.152 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.154 

0.154 

0.154 

0.155 

0.155 

0.155 



75 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.153 

0.154 

0.154 

0.155 

0.155 

0.156 




80 

0.153 

0.153 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.153 

0.153 

0.153 

0.154 

0.154 

0.155 





85 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.153 

0.153 

0.152 

0.152 

0.152 

0.152 






90 

0.154 

0.154 

0.154 

0.154 

0.154 

0.154 

0.153 

0.152 

0.151 

0.150 

0.149 

0.148 

0.147 








Table 6.— F,(a,Q 


20 

0 

5 

10 

15 

20 

25 

30 

35 

40 

0 

0.141 

0.141 

0.141 

0.141 

0.142 

0.142 

0.142 

0.143 

0.143 

5 

0.141 

0.141 

0.141 

0.141 

0.142 

0.142 

0.142 

0.143 

0.143 

10 

0.141 

0.141 

0.141 

0.142 

0.142 

0.142 

0.143 

0.143 

0.144 

15 

0.141 

0.141 

0.142 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

20 

0.142 

0.142 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

0.144 

25 

0.142 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

0.144 

0.145 

30 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

0.144 

0.145 

0.145 

35 

0.143 

0.143 

0.143 

0.143 

0.144 

0.144 

0.145 

0.145 

0.145 

40 

0.143 

0.143 

0.144 

0.144 

0.144 

0.145 

0.145 

0.145 

0.146 

45 

0.144 

0.144 

0.144 

0.144 

0.145 

0.145 

0.145 

0.146 

0.146 

50 

0.144 

0.144 

0.144 

0.145 

0.145 

0.146 

0.146 

0.146 

0.147 

55 

0.144 

0.145 

0.145 

0.145 

0.146 

0.146 

0.146 

0.147 

0.147 

60 

0.145 

0.145 

0.145 

0.146 

0.146 

0.147 

0.147 

0.147 

0.148 

65 

0.145 

0.145 

0.146 

0.146 

0.147 

0.147 

0.147 

0.148 

0.148 

70 

0.146 

0.146 

0.146 

0.147 

0.147 

0.147 

0.148 

0.148 

0.148 

75 

0.146 

0.146 

0.146 

0.147 

0.147 

0.147 

0.147 

0.147 

0.148 

80 

0.145 

0.145 

0.146 

0.146 

0.146 

0.147 

0.147 

0.146 

0.146 

85 

0.144 

0.144 

0.145 

0.145 

0.145 

0.145 

0.145 

0.144 

0.144 

90 

0.142 

0.142 

0.143 

0.143 

0.143 

0.142 

0.141 

0.140 

0.139 


i, 02 ) for a = 20. 


45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

0.144 

0.144 

0.144 

0.145 

0.145 

0.146 

0.146 

0.145 

0.144 

0.142 

0.144 

0.144 

0.145 

0.145 

0.145 

0.146 

0.146 

0.145 

0.144 

0.142 

0.144 

0.144 

0.145 

0.145 

0.146 

0.146 

0.146 

0.146 

0.145 

0.143 

0.144 

0.145 

0.145 

0.146 

0.146 

0.147 

0.147 

0.146 

0.145 

0.143 

0.145 

0.145 

0.146 

0.146 

0.147 

0.147 

0.147 

0.146 

0.145 

0.143 

0.145 

0.146 

0.146 

0.147 

0.147 

0.147 

0.147 

0.147 

0.145 

0.142 

0.145 

0.146 

0.146 

0.147 

0.147 

0.148 

0.147 

0.147 

0.145 

0.141 

0.146 

0.146 

0.147 

0.147 

0.148 

0.148 

0.147 

0.146 

0.144 

0.140 

0.146 

0.147 

0.147 

0.148 

0.148 

0.148 

0.148 

0.146 

0.144 

0.139 

0.147 

0.147 

0.148 

0.148 

0.148 

0.148 

0.148 

0.146 

0.143 

0.137 

0.147 

0.148 

0.148 

0.148 

0.149 

0.149 

0.148 

0.146 

0.143 

0.136 

0.148 

0.148 

0.148 

0.149 

0.149 

0.149 

0.149 

0.147 

0.142 


0.148 

0.148 

0.149 

0.149 

0.150 

0.150 

0.149 

0.147 



0.148 

0.149 

0.149 

0.150 

0.150 

0.151 

0.151 




0.148 

0.149 

0.149 

0.150 

0.151 

0.152 





0.148 

0.148 

0.149 

0.149 

0.151 






0.146 

0.146 

0.147 

0.147 







0.143 

0.143 

0.142 








0.137 

0.136 










Table 7. — 7 7 ,{a,0|,0 2 ) for a = 25. 


NASA/CR— 2005-213799 


20 



25 

0 

5 


15 


25 





50 

55 


65 


75 

80 

85 

90 

0 

0.136 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.139 

0.139 

0.139 

0.140 

0.140 

0.140 

0.140 

0.139 

0.138 

0.136 

0.133 

0.129 

5 

0.137 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.140 

0.140 

0.140 

0.140 

0.138 

0.137 

0.134 

0.129 

10 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.140 

0.141 

0.141 

0.141 

0.140 

0.139 

0.137 

0.134 

0.130 

15 

0.137 

0.137 

0.138 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.141 

0.141 

0.141 

0.141 

0.140 

0.138 

0.135 

0.130 

20 

0.138 

0.138 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.141 

0.142 

0.142 

0.142 

0.141 

0.140 

0.138 

0.135 

0.130 

25 

0.138 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.142 

0.142 

0.142 

0.142 

0.141 

0.139 

0.135 

0.129 

30 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.142 

0.143 

0.143 

0.143 

0.142 

0.141 

0.139 

0.135 

0.128 

35 

0.139 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

0.143 

0.143 

0.141 

0.139 

0.134 

0.127 

40 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.143 

0.144 

0.144 

0.144 

0.143 

0.142 

0.139 

0.133 

0.125 

45 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.144 

0.144 

0.144 

0.145 

0.145 

0.144 

0.142 

0.139 

0.133 


50 

0.140 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.144 

0.144 

0.145 

0.145 

0.145 

0.145 

0.143 

0.139 



55 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.144 

0.144 

0.145 

0.146 

0.146 

0.146 

0.146 

0.144 




60 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.144 

0.144 

0.145 

0.145 

0.146 

0.147 

0.147 

0.147 





65 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.144 

0.145 

0.145 

0.146 

0.147 

0.148 






70 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.143 

0.143 

0.144 

0.145 

0.146 

0.147 







75 

0.138 

0.138 

0.139 

0.140 

0.140 

0.141 

0.141 

0.141 

0.142 

0.142 

0.143 

0.144 








80 

0.136 

0.137 

0.137 

0.138 

0.138 

0.139 

0.139 

0.139 

0.139 

0.139 

0.139 









85 

0.133 

0.134 

0.134 

0.135 

0.135 

0.135 

0.135 

0.134 

0.133 

0.133 










90 

0.129 

0.129 

0.130 

0.130 

0.130 

0.129 

0.128 

0.127 

0.125 












30 

0 

5 

10 

15 

20 

25 

Table 8 . — for a = 

30 35 40 45 50 

30. 

55 

60 

65 

70 

75 

80 

85 

90 

0 

0.133 

0.133 

0.133 

0.134 

0.134 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.135 

0.134 

0.133 

0.130 

0.127 

0.122 

0.116 

5 

0.133 

0.133 

0.134 

0.134 

0.135 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.136 

0.135 

0.133 

0.131 

0.128 

0.123 

0.117 

10 

0.133 

0.134 

0.134 

0.134 

0.135 

0.135 

0.136 

0.136 

0.137 

0.137 

0.137 

0.137 

0.137 

0.136 

0.134 

0.132 

0.129 

0.124 

0.118 

15 

0.134 

0.134 

0.134 

0.135 

0.135 

0.136 

0.136 

0.137 

0.137 

0.138 

0.138 

0.138 

0.137 

0.137 

0.135 

0.133 

0.130 

0.125 

0.119 

20 

0.134 

0.135 

0.135 

0.135 

0.136 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.138 

0.138 

0.138 

0.136 

0.134 

0.130 

0.125 

0.119 

25 

0.135 

0.135 

0.135 

0.136 

0.137 

0.137 

0.138 

0.138 

0.139 

0.139 

0.139 

0.139 

0.139 

0.138 

0.137 

0.134 

0.130 

0.125 

0.118 

30 

0.135 

0.135 

0.136 

0.136 

0.137 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.140 

0.140 

0.139 

0.137 

0.135 

0.131 

0.125 

0.117 

35 

0.136 

0.136 

0.136 

0.137 

0.137 

0.138 

0.139 

0.139 

0.140 

0.140 

0.140 

0.141 

0.140 

0.140 

0.138 

0.135 

0.130 

0.124 


40 

0.136 

0.136 

0.137 

0.137 

0.138 

0.139 

0.139 

0.140 

0.140 

0.141 

0.141 

0.141 

0.141 

0.140 

0.139 

0.136 

0.131 



45 

0.136 

0.136 

0.137 

0.138 

0.138 

0.139 

0.140 

0.140 

0.141 

0.141 

0.142 

0.142 

0.142 

0.141 

0.140 

0.136 




50 

0.136 

0.136 

0.137 

0.138 

0.138 

0.139 

0.140 

0.140 

0.141 

0.142 

0.142 

0.143 

0.143 

0.142 

0.141 





55 

0.136 

0.136 

0.137 

0.138 

0.138 

0.139 

0.140 

0.141 

0.141 

0.142 

0.143 

0.143 

0.144 

0.144 






60 

0.135 

0.136 

0.137 

0.137 

0.138 

0.139 

0.140 

0.140 

0.141 

0.142 

0.143 

0.144 

0.145 







65 

0.134 

0.135 

0.136 

0.137 

0.138 

0.138 

0.139 

0.140 

0.140 

0.141 

0.142 

0.144 








70 

0.133 

0.133 

0.134 

0.135 

0.136 

0.137 

0.137 

0.138 

0.139 

0.140 

0.141 









75 

0.130 

0.131 

0.132 

0.133 

0.134 

0.134 

0.135 

0.135 

0.136 

0.136 










80 

0.127 

0.128 

0.129 

0.130 

0.130 

0.130 

0.131 

0.130 

0.131 











85 

0.122 

0.123 

0.124 

0.125 

0.125 

0.125 

0.125 

0.124 












90 

0.116 

0.117 

0.118 

0.119 

0.119 

0.118 

0.117 














35 

0 

5 

10 

15 

20 

25 

Table 9. — -i 7 J {a,0i,0 2 ) for a = 

30 35 40 45 50 

35. 

55 

60 

65 

70 

75 

80 

85 

90 

0 

0.130 

0.131 

0.131 

0.132 

0.132 

0.132 

0.133 

0.133 

0.133 

0.133 

0.133 

0.132 

0.131 

0.129 

0.127 

0.123 

0.118 

0.1 

113 

0.108 

5 

0.131 

0.131 

0.131 

0.132 

0.132 

0.133 

0.133 

0.133 

0.134 

0.134 

0.133 

0.133 

0.132 

0.130 

0.127 

0.124 

0.120 

0.1 

115 

0.111 

10 

0.131 

0.131 

0.132 

0.132 

0.133 

0.133 

0.134 

0.134 

0.134 

0.134 

0.134 

0.134 

0.133 

0.131 

0.129 

0.126 

0.121 

0.1 

117 

0.113 

15 

0.132 

0.132 

0.132 

0.133 

0.133 

0.134 

0.134 

0.135 

0.135 

0.135 

0.135 

0.135 

0.134 

0.132 

0.130 

0.127 

0.123 

0.1 

118 

0.114 

20 

0.132 

0.132 

0.133 

0.133 

0.134 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.135 

0.134 

0.131 

0.128 

0.123 

0.1 

118 

0.113 

25 

0.132 

0.133 

0.133 

0.134 

0.135 

0.135 

0.136 

0.136 

0.137 

0.137 

0.137 

0.137 

0.136 

0.135 

0.132 

0.128 

0.123 

0.1 

118 


30 

0.133 

0.133 

0.134 

0.134 

0.135 

0.136 

0.137 

0.137 

0.138 

0.138 

0.138 

0.138 

0.137 

0.136 

0.133 

0.129 

0.123 




35 

0.133 

0.133 

0.134 

0.135 

0.136 

0.136 

0.137 

0.138 

0.138 

0.139 

0.139 

0.139 

0.138 

0.136 

0.134 

0.129 





40 

0.133 

0.134 

0.134 

0.135 

0.136 

0.137 

0.138 

0.138 

0.139 

0.139 

0.140 

0.140 

0.139 

0.138 

0.135 






45 

0.133 

0.134 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.139 

0.140 

0.141 

0.141 

0.140 

0.139 







50 

0.133 

0.133 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.141 

0.141 

0.142 

0.141 








55 

0.132 

0.133 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.141 

0.142 

0.142 









60 

0.131 

0.132 

0.133 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.141 










65 

0.129 

0.130 

0.131 

0.132 

0.134 

0.135 

0.136 

0.136 

0.138 

0.139 











70 

0.127 

0.127 

0.129 

0.130 

0.131 

0.132 

0.133 

0.134 

0.135 












75 

0.123 

0.124 

0.126 

0.127 

0.128 

0.128 

0.129 

0.129 













80 

0.118 

0.120 

0.121 

0.123 

0.123 

0.123 

0.123 














85 

0.113 

0.115 

0.117 

0.118 

0.118 

0.118 















90 

0.108 

0.111 

0.113 

0.114 

0.113 

















Table 10. — -F/a, 61 , 62 ) for a = 40. 
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I 25 30 35 40 45 50 55 60 65 70 75 80 


0 

0.129 

0.129 

0.129 

0.130 

0.131 

0.131 

0.131 

0.132 

0.131 

0.131 

0.130 

0.129 

0.127 

0.124 

0.121 

0.117 

0.113 

5 

0.129 

0.129 

0.130 

0.130 

0.131 

0.131 

0.132 

0.132 

0.132 

0.132 

0.131 

0.130 

0.128 

0.126 

0.122 

0.119 

0.115 

10 

0.129 

0.130 

0.130 

0.131 

0.132 

0.132 

0.133 

0.133 

0.133 

0.133 

0.132 

0.131 

0.130 

0.127 

0.124 

0.121 

0.118 

15 

0.130 

0.130 

0.131 

0.132 

0.132 

0.133 

0.134 

0.134 

0.134 

0.134 

0.134 

0.133 

0.131 

0.129 

0.126 

0.123 

0.120 

20 

0.131 

0.131 

0.132 

0.132 

0.133 

0.134 

0.134 

0.135 

0.135 

0.135 

0.135 

0.134 

0.133 

0.130 

0.127 

0.124 

0.120 

25 

0.131 

0.131 

0.132 

0.133 

0.134 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.135 

0.134 

0.132 

0.128 

0.124 


30 

0.131 

0.132 

0.133 

0.134 

0.134 

0.135 

0.136 

0.137 

0.137 

0.137 

0.137 

0.137 

0.135 

0.133 

0.129 



35 

0.132 

0.132 

0.133 

0.134 

0.135 

0.136 

0.137 

0.137 

0.138 

0.138 

0.138 

0.138 

0.136 

0.134 




40 

0.131 

0.132 

0.133 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.139 

0.139 

0.139 

0.138 





45 

0.131 

0.132 

0.133 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.140 

0.140 






50 

0.130 

0.131 

0.132 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.141 







55 

0.129 

0.130 

0.131 

0.133 

0.134 

0.135 

0.137 

0.138 

0.139 

0.140 








60 

0.127 

0.128 

0.130 

0.131 

0.133 

0.134 

0.135 

0.136 

0.138 









65 

0.124 

0.126 

0.127 

0.129 

0.130 

0.132 

0.133 

0.134 










70 

0.121 

0.122 

0.124 

0.126 

0.127 

0.128 

0.129 











75 

0.117 

0.119 

0.121 

0.123 

0.124 

0.124 












80 

0.113 

0.115 

0.118 

0.120 

0.120 













85 

0.111 

0.114 

0.118 

0.119 














90 

0.112 

0.119 

0.123 
















0.111 

0.114 

0.118 

0.119 


Table 11. — ■F I {a,0i,0 2 ) for a = 45. 


45 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

0 

0.128 

0.128 

0.129 

0.129 

0.130 

0.130 

0.131 

0.131 

0.130 

0.130 

0.128 

0.126 

0.124 

0.121 

0.117 

0.115 

0.116 

5 

0.128 

0.128 

0.129 

0.130 

0.130 

0.131 

0.131 

0.131 

0.131 

0.130 

0.129 

0.127 

0.125 

0.123 

0.120 

0.119 

0.123 

10 

0.129 

0.129 

0.130 

0.131 

0.131 

0.132 

0.132 

0.132 

0.132 

0.132 

0.131 

0.130 

0.128 

0.125 

0.123 

0.123 

0.129 

15 

0.129 

0.130 

0.131 

0.131 

0.132 

0.133 

0.133 

0.134 

0.134 

0.134 

0.133 

0.132 

0.130 

0.127 

0.126 

0.125 


20 

0.130 

0.130 

0.131 

0.132 

0.133 

0.134 

0.135 

0.135 

0.135 

0.135 

0.134 

0.133 

0.131 

0.129 

0.127 



25 

0.130 

0.131 

0.132 

0.133 

0.134 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.135 

0.133 

0.130 




30 

0.131 

0.131 

0.132 

0.133 

0.135 

0.136 

0.136 

0.137 

0.137 

0.138 

0.137 

0.136 

0.134 





35 

0.131 

0.131 

0.132 

0.134 

0.135 

0.136 

0.137 

0.138 

0.138 

0.139 

0.138 

0.137 






40 

0.130 

0.131 

0.132 

0.134 

0.135 

0.136 

0.137 

0.138 

0.139 

0.140 

0.140 







45 

0.130 

0.130 

0.132 

0.134 

0.135 

0.136 

0.138 

0.139 

0.140 

0.141 








50 

0.128 

0.129 

0.131 

0.133 

0.134 

0.136 

0.137 

0.138 

0.140 









55 

0.126 

0.127 

0.130 

0.132 

0.133 

0.135 

0.136 

0.137 










60 

0.124 

0.125 

0.128 

0.130 

0.131 

0.133 

0.134 











65 

0.121 

0.123 

0.125 

0.127 

0.129 

0.130 












70 

0.117 

0.120 

0.123 

0.126 

0.127 













75 

0.115 

0.119 

0.123 

0.125 














80 

0.116 

0.123 

0.129 















85 

0.124 

0.140 
















90 

0.166 


















Table 12. — i 7 J (a,0i,0 2 ) for a = 50. 


50 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

0 

0.127 

0.128 

0.129 

0.129 

0.130 

0.131 

0.131 

0.130 

0.130 

0.128 

0.127 

0.124 

0.122 

0.120 

0.120 

5 

0.128 

0.128 

0.129 

0.130 

0.131 

0.131 

0.131 

0.131 

0.131 

0.130 

0.128 

0.126 

0.125 

0.124 

0.128 

10 

0.129 

0.129 

0.130 

0.131 

0.132 

0.132 

0.133 

0.133 

0.133 

0.132 

0.131 

0.129 

0.128 

0.129 

0.133 

15 

0.129 

0.130 

0.131 

0.132 

0.133 

0.134 

0.134 

0.135 

0.135 

0.134 

0.133 

0.132 

0.131 

0.131 


20 

0.130 

0.131 

0.132 

0.133 

0.134 

0.135 

0.136 

0.136 

0.136 

0.136 

0.135 

0.134 

0.132 



25 

0.131 

0.131 

0.132 

0.134 

0.135 

0.136 

0.137 

0.138 

0.138 

0.137 

0.137 

0.135 




30 

0.131 

0.131 

0.133 

0.134 

0.136 

0.137 

0.138 

0.139 

0.139 

0.139 

0.138 





35 

0.130 

0.131 

0.133 

0.135 

0.136 

0.138 

0.139 

0.140 

0.140 

0.140 






40 

0.130 

0.131 

0.133 

0.135 

0.136 

0.138 

0.139 

0.140 

0.141 







45 

0.128 

0.130 

0.132 

0.134 

0.136 

0.137 

0.139 

0.140 








50 

0.127 

0.128 

0.131 

0.133 

0.135 

0.137 

0.138 









55 

0.124 

0.126 

0.129 

0.132 

0.134 

0.135 










60 

0.122 

0.125 

0.128 

0.131 

0.132 











65 

0.120 

0.124 

0.129 

0.131 












70 

0.120 

0.128 

0.133 














75 0.126 0.143 

80 0.166 


75 80 

0.126 0.166 
0.143 


Table 13. — ■i r I {a,0i,0 2 ) for a = 55. 


0.112 

0.119 

0.123 


90 

0.166 


NASA/CR— 2005-213799 


22 



0 

0.128 

0.128 

0.129 

0.130 

0.131 

0.131 

0.131 

0.131 

0.129 

0.128 

0.126 

0.124 

0.124 

5 

0.128 

0.129 

0.130 

0.131 

0.132 

0.132 

0.132 

0.132 

0.131 

0.130 

0.129 

0.129 

0.132 

10 

0.129 

0.130 

0.131 

0.133 

0.133 

0.134 

0.134 

0.134 

0.134 

0.133 

0.133 

0.134 

0.138 

15 

0.130 

0.131 

0.133 

0.134 

0.135 

0.136 

0.136 

0.137 

0.136 

0.136 

0.136 

0.136 


20 

0.131 

0.132 

0.133 

0.135 

0.136 

0.137 

0.138 

0.139 

0.139 

0.138 

0.138 



25 

0.131 

0.132 

0.134 

0.136 

0.137 

0.139 

0.140 

0.140 

0.140 

0.140 




30 

0.131 

0.132 

0.134 

0.136 

0.138 

0.140 

0.141 

0.141 

0.142 





35 

0.131 

0.132 

0.134 

0.137 

0.139 

0.140 

0.141 

0.142 






40 

0.129 

0.131 

0.134 

0.136 

0.139 

0.140 

0.142 







45 

0.128 

0.130 

0.133 

0.136 

0.138 

0.140 








50 

0.126 

0.129 

0.133 

0.136 

0.138 









55 

0.124 

0.129 

0.134 

0.136 










60 

0.124 

0.132 

0.138 











65 

0.129 

0.146 












70 

0.166 














60 

0 

5 

10 

Table 14.- 

15 20 

-F/(a, 0 i, 0 2 ) for a 

25 30 35 

= 60. 
40 

45 

50 

55 

0 

0.129 

0.130 

0.131 

0.132 

0.133 

0.133 

0.132 

0.131 

0.130 

0.128 

0.128 

0.131 

5 

0.130 

0.131 

0.132 

0.133 

0.134 

0.134 

0.134 

0.134 

0.133 

0.133 

0.136 

0.149 

10 

0.131 

0.132 

0.134 

0.135 

0.136 

0.137 

0.137 

0.137 

0.137 

0.138 

0.142 


15 

0.132 

0.133 

0.135 

0.137 

0.138 

0.139 

0.140 

0.140 

0.140 

0.141 



20 

0.133 

0.134 

0.136 

0.138 

0.140 

0.141 

0.142 

0.142 

0.143 




25 

0.133 

0.134 

0.137 

0.139 

0.141 

0.142 

0.143 

0.144 





30 

0.132 

0.134 

0.137 

0.140 

0.142 

0.143 

0.145 






35 

0.131 

0.134 

0.137 

0.140 

0.142 

0.144 







40 

0.130 

0.133 

0.137 

0.140 

0.143 








45 

0.128 

0.133 

0.138 

0.141 









50 

0.128 

0.136 

0.142 










55 

0.131 

0.149 











60 

0.166 













Table 15. — F,(a,0 1 ,0 2 ) for a = 65. 


65 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

0 

0.131 

0.132 

0.134 

0.135 

0.135 

0.135 

0.134 

0.132 

0.132 

0.134 

5 

0.132 

0.133 

0.135 

0.136 

0.137 

0.137 

0.137 

0.138 

0.140 

0.152 

10 

0.134 

0.135 

0.137 

0.139 

0.140 

0.141 

0.142 

0.143 

0.147 


15 

0.135 

0.136 

0.139 

0.141 

0.143 

0.144 

0.145 

0.146 



20 

0.135 

0.137 

0.140 

0.143 

0.145 

0.146 

0.147 




25 

0.135 

0.137 

0.141 

0.144 

0.146 

0.148 





30 

0.134 

0.137 

0.142 

0.145 

0.147 






35 

0.132 

0.138 

0.143 

0.146 







40 

0.132 

0.140 

0.147 








45 

0.134 

0.152 









50 

0.166 











70 

0 

5 

10 

Table 16.- 

15 20 

-F / (a,0 1 ,0 2 ) for a = 70. 

25 30 35 40 

0 

0.133 

0.135 

0.137 

0.138 

0.138 

0.137 

0.136 

0.136 0.166 

5 

0.135 

0.137 

0.139 

0.141 

0.142 

0.142 

0.145 

0.155 

10 

0.137 

0.139 

0.142 

0.145 

0.146 

0.148 

0.151 


15 

0.138 

0.141 

0.145 

0.147 

0.149 

0.151 



20 

0.138 

0.142 

0.146 

0.149 

0.152 




25 

0.137 

0.142 

0.148 

0.151 





30 

0.136 

0.145 

0.151 






35 

0.136 

0.155 







40 

0.166 









Table 17. — F, {a, 61 , 62 ) for a = 75. 
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75 

0 

5 

10 

15 

20 

25 

30 

0 

0.136 

0.138 

0.141 

0.141 

0.140 

0.139 

0.166 

5 

0.138 

0.142 

0.145 

0.147 

0.149 

0.158 


10 

0.141 

0.145 

0.149 

0.153 

0.156 



15 

0.141 

0.147 

0.153 

0.156 




20 

0.140 

0.149 

0.156 





25 

0.139 

0.158 






30 

0.166 








Table 18. — F, (a, 61 , 62 ) for a = 80. 


80 

0 

5 

10 

15 

20 

0 

0.139 

0.143 

0.145 

0.143 

0.166 

5 

0.143 

0.149 

0.154 

0.161 


10 

0.145 

0.154 

0.161 



15 

0.143 

0.161 




20 

0.166 






Table 19. — F t (a,0i,0 2 ) for a = 85. 


85 

0 

5 

10 

0 

0.142 

0.148 

0.166 

5 

0.148 

0.165 


10 

0.166 




Table 20. — -F,(a,0|,0 2 ) for a = 90. 


90 

0 

0 

0.167 
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(or more) contact angles of the system. A numerical solution is employed to verify the new problem formulation. The benchmarked 
computations support the use of the existing theoretical approach to geometries with variable wettability. Simple experiments to 
confirm the theoretical findings are recommended. Favorable agreement between such experiments and the present theory may argue 
well for the extension of the analytic results to predict fluid performance in future large length scale capillary fluid systems for 
spacecraft as well as for small scale capillary systems on Earth. 
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